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INTEGRAL PRODUCTS AND PROBABILITY. 

By P. J. Daniell. 

1. Introduction. — In many problems arising in statistical biology and 
statistical economics time enters as an indispensable factor. It is the chief 
aim of this paper to provide a form of analysis suitable for such problems, 
and this will be a theory of probability using time as an auxiliary variable. 
The term " Dynamic Probability " may be used to avoid confusion with 
Gibbs' Statistical Mechanics, which is a theory of interactions between 
fast-moving molecules or between vibratory systems, whereas our theory 
is that of time- variations in position of groups moving more or less erratic- 
ally. That is to say, the basis is change of position rather than of velocity, 
and the tendency to disperse is regarded as inherent, not as the effect of 
" collisions." 

For the sake of simplicity we consider only motion in a single dimension. 
The variable x may be a measure of actual geometrical position, or it may 
measure some factor of the relative environment, for example, temperature, 
intensity of light, or it may refer to a quantity of goods. 

The first step in the analysis is a search for some standard formula on 
which may be built a more complex and general theory. It is found that, 
if certain natural assumptions are made, a functional equation is satisfied, 
which is expressed in terms of a Stieltjes integral product. This form of 
product may be compared with Volterra's integral composition, and on it a 
similar algebra may be built. It is shown that there exists an idem-factor 
function for this algebra whereas such a function does not exist for the 
Volterra composition.* The Stieltjes integral product itself forms a 
secondary nucleus for our paper; its investigation will be found in para- 
graph 3. 

2. Fundamental Equation. — Denote an increase (possibly negative) in x 
by y and an increase in time t by u. The motion in the interval u will 
depend partly on the previous history H of the member considered. Denote 

by 
P(yu «i, Ho) 

*V. Volterra, "Lemons sur les fonctions de lignes." Paris (1913). G. C. Evans, 
Cambridge Colloquium, New York (1916), p. 117. By introducing an auxiliary algebraic 
symbol j, Evans modifies the Volterra algebra so as to contain an idem-factor but this 
factor, 1 + jO, is independent of the variables. For a resume of the Stieltjes integral and 
references the reader is referred to T. H. Hildebrandt, Bulletin of the American Mathematical 
Society (1917), 24, p. 117 and p. 177. 
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144 Daniell: Integral Products and Probabilities. 

the probability that, in an interval u\, the increase in x does not exceed 
yi (i.e., y ^ yi) for a certain history H up to time to. Similarly we may 
denote by 

P(yi, u 2 , Hi) 

the probability that, in a succeeding interval « 2 , the increase in x does not 
exceed y 2 . Hi will include both H and the motion during the interval «i. 

This is the critical point at which the dynamics of material and living 
objects diverge. Theoretical mechanics is expressed in terms of velocities 
and accelerations which are instantaneous time-derivatives. In such a 
theory, however small we may choose the interval u\, the second probability 
P(y, ui, Hi) depends essentially on the motion during ui rather than on H . 
If only Mi is chosen sufficiently small the contrary is true of living objects. 
It is to be understood that we are not concerned with the ultimate 
hypotheses of materialism and vitalism; we regard each object as a whole 
and do not pretend to delve into the dynamics of its finest constituents 
which may or may not be of the material type. Our point will be made 
clear by an illustration. Consider a salmon ascending a stream and 
suddenly brought to rest by a glass plate, inserted in the stream for a brief 
instant. After the plate is removed the salmon will resume its motion 
almost as if there had been no interruption. It will be influenced more by 
previous history or habit than by the interruption, provided the latter is 
short. Compare this case with that of a cylinder rolling along a horizontal 
plane. Even after the very shortest interruption the cylinder will remain 
at rest. It is more influenced by the immediately preceding history. An 
example which is an apparent exception is even more illuminating. In a 
speculative stock market the movement of a particular stock depends 
partly on the preceding movement. But even here the effect is influenced 
by the events of some hours or minutes rather than by those of the last 
hundredth of a second, let us say. 

We shall assume that if the interval ui is sufficiently small, the probability 
function P depends on Ho only, so that 

P(y 2 , u 2 , Hi) = P{y 2 , u^., H ). 

If we confine our attention to a homogeneous group of members having 
practically the same history H , up to the time to, then we may omit the 
symbol H and define P(y, u) as the probability that the increase in x 
does not exceed y in a short interval of time u not far removed from a 
fixed point of time t . 

This is the fundamental assumption on which our analysis is based; 
and it is this which sharply distinguishes dynamic probability from the 
classical theories of mechanics, statistical or otherwise. After an investiga- 
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tion of P{y, u) in which u is assumed to be small it is possible to build up a 
more general theory suitable for longer intervals, somewhat as organic 
chemistry is built on the atomic laws of simpler reactions. 

By its definition P(y, u) is a limited non-decreasing function of y ap- 
proaching as y approaches — oo , and 1 as y approaches + <x> . It is also 
" continuous on the right," that is 

P(y, u) = lim P(y + e, u) 

e=0, e>0 

= P(m+0,m). 
On the other hand, 

P(y — 0, u) — lim P(y — e, u) 

is the probability that the increase in x is strictly less than y, and it is not 
necessarily equal to P(y, u). 

The increase y in the interval u x + m 2 is the sum of an increase z during Mi 
and an increase y — z during m 2 . Keeping y fixed the function P(y — z, u 2 ) 
is a limited non-increasing function of z continuous on the left with the 
values 1 at — oo , at + <» , and is continuous in each of a countable set of. 
intervals which are open below and closed above. 

Let cd be one of these intervals; then the probability that, in the time 
Mi, the increase z in x satisfies c < z ^d will be 

AP(z, u{) = P(d, tii) — P(c, ui). 

Given that c < z ^ d, the probability that the increase in the time u\ -f n? 
does not exceed y will lie between 

P(y — d, u 2 ), P(y — c, u 2 ). 

Then the combined probability will lie between 

P(y - d, u 2 )AP(z, Mi), P(y — c, u 2 )AP(z, Mi). 

The law of composite probability is applicable here because we assumed 
the motion during u 2 independent of that in u\. 

Divide cd into smaller sub-intervals and proceed to the limit as their 
maximum length approaches 0. The combined probability that c < z ^ d 
and that during Mi + m 2 the increase does not exceed y will be 



r 



P(y — z, u 2 )d z P(z, Mi). 



Add up the corresponding expressions for all the intervals cd, and then the 
fundamental equation is obtained, namely 

2(1) P(y, Mi + «,) = f " P(y - z, u 2 )d z P(z, Mi). 
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Without changes in the processes of reasoning we may regard u?, as the 
earlier, Mi the later interval of time and so obtain 

2(2) P(y, Mi + u 2 ) = f P(y - z, wiKPfe u*). 

J— 00 

3. Stieltjes Integral Products. — The function a(x) will be said to be a 
regular function of limited variation in x if it is of limited variation 
(— oo to + oo ), if a(— oo ) = and if, for all values of x, 

a(x) = \[a{x + 0) + a(x - 0)]. 

We shall require an important lemma. 

3(1). If a(s, t) is a regular function of limited variation in s, and measur- 
able in the sense of Borel in t, and if the variation be limited uniformly with 
respect to t, and if g(t) is of limited variation and f(s) is bounded and 
measurable,* then 

r"f(«)d..r"a(s, t)gd(t) = f + " f" f + "/(*)<W*, f)~\dg(t). 

J—at> */— oo *J— oo L*' — oo J 

Let/(s) = 1, (c ^ s ^ d) = otherwise. 

Cf(»)d. r +C °«(5, t)dg(t) = f *a(d+0, t)dg(t) - f "a(c - 0, 0<^(0 

X+oo r /»+« ' j 

LJ-oo /(*)*"(*'*) J ^O- 

If /(s) is a step-function (constant over each of a finite set of sub-intervals) 
it is a linear combination of functions of the above type and therefore the 
equality will be satisfied. But any measurable function f(s) can be obtained 
from such step-functions by a finite succession of linear combinations and 
limiting processes using in all cases bounded sets of functions since f(s) 
is bounded. Because a(s, t) is of uniformly limited variation a number 8 
exists such that the total variation of a(s, t) does not exceed S. If 6 is the 

X+oo 
a(s, t)dg(t) does not exceed SG. 

Also if |/(5)| <M, 

f + °7(*Ka(s, t) < MS. 

I t/ — 00 i 

Therefore in each of the required limiting processes the several integrands 

do not exceed summable functions in absolute value, and at each step the 

equality is preserved, f 

* In this paper the attribute 'measurable' will be taken in the sense of Borel. 
t P. J. Daniell, Annals of Mathematics, Vol. 19 (1918), p. 290. The theorem used here 
is used frequently in the present paper. 
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3(2). Let ai(x, y), a 2 (x, y) be regular functions of limited variation in x, 
and measurable in y, and let their total variations be limited uniformly 
with respect to y; then we define their 8 — V composition (Stieltjes- 
Volterra integral product) as 

ai(x, s)d„a 2 (s, y). Def. 

00 

Evidently this " product " satisfies the distributive laws 

(«1 + <X 2 )-a 3 = «l-«3 + Ol2-Ot3, 

ai-(«2 + « 3 ) = «i-a 2 + ai-a 3 . 

In theorem 3(1) replace /(s) by ai(x, s), a(s, t) by a 2 {s, t) and g(t) by a 3 (t, y). 
The theorem becomes 

I ai(x, s)d 3 J a 2 {s, t)d t a 3 (t, y) 

I ai(x, s)d s a 2 (s, t) d t a 3 (t, y), 

which is, symbolically, the associative law 

ai-{a 2 -<Xi) = (ai-a 2 )-a z . 
Consider the function t?(x, y) = 0, x < y, 



— 2> % — y> 



= 1, x > y. 

This satisfies the conditions imposed on «i, a 2 , and 

/•+00 

«•»?(«, 2/) = I «(a;, sRt?(s, 1/) = a(x, y), 
J— K 

77(0:, s)d s a{s, y) 

00 

= a(x - 0, 2/) - a(- 00, 7/) + |[a(x + 0, ?/) - a(a; - 0, 2/)] 

= «(«> */), 

since a is a regular function of limited variation in x. We see that 77 plays 
the part of an idem-factor in either order, 

a- 77 = 7] -a = a. 

Now let 77' be another idem-factor satisfying the conditions in 3(2) such that 
for every a, 77' -a = a. Then by what has been proved 77' • 77 = 77' and by 
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hypothesis v'-v = r\ so that 77' = n. Similarly it can be shown that there 
is no other idem-fector r\" satisfying the conditions in 3(2) and such that 
a • 77" = a for every a. 

Thus 77 is the only idem-factor in the class of functions a considered. 

If it happens that a\-a 2 = a 2 -a\ we may say that the functions are 
S-permutable to distinguish the property from that of being permutable 
according to Volterra. 

If a u a 2 are differentiable with respect to their first variables, for example 



then 



fi(s, t)ds, a 2 (s, t) = I f 2 (s, t)ds, 

00 %) — 00 



X' * * 
00 



t)ds, 



in the Volterra notation, and this shows the relation between the two types 
of combination. 

A development of this algebra particularly in connection with integral 
equations would be interesting but we shall now confine our attention to 
functions of the difference x — y, 

a(x, y) = a(x - y). 

Then 

X+oo /»+» 

a x (x — s)d s a 2 (s — y) = I ai(x — y — z)d z a 2 (z), 
00 *y—ao 

so that a\-a 2 is also a function of the difference x — y. 

3(3). To use symbols more in keeping with our special problem, if 
<xi(y)> <x 2 (y) are regular functions of limited variation (— 00 to + 00), then 
we define their S — V product as 



«i 



a x (y - z)da 2 (z) Def. 

00 



The corresponding idem-factor will be 

v(y) = 0, y < 0, 

= h y = o, 
= 1, y > 0. 

3(4). Theorem. The product as defined in 3(3) is S-per •mutable. For 
since the functions are regular, by a theorem on integration by parts* 

* P. J. Daniell, Transactions of the American Mathematical Society, Vol. 19 (1918), 
p. 362. 
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ai(y — z)da 2 (z) 

00 

/>+00 
= [«i(2/ — 2)«2(z)]±2 — / (x 2 (z)d z ai(y — z). 

«/ — X 

Butai(2/— co),a 2 (— <») are each while a 2 (+ °o),ai(y+ °o) are finite. 
Again if w — z = v, 

kzfxiiy — z) = «i(w — z 2 ) — ai(t/ — 2i), zi < z 2 , 

= ai(r 2 ) — ai(»i), v x > v 2 , 

= — A„ai(»). 

Combining these results, 

oci-a 2 (y) = I a 2 (y — v)doti(v) = a 2 -ai(y). 

This proves the theorem. 

In the fundamental equation 2(1) P(y, u) is not regular, but we can 
replace it by ft,(y) = \[P{y + 0, u) + P(y - 0, «)]. 

Equation 2(1) can now be written 

X+00 
P(w - z, u 2 )dp ui (z), 
00 

or what is the same equation, 

P(y + 0, mi + «») = f P(y + - z, u 2 )d(3 Ul (z). 

J — 00 

Again 

X+00 
P(y — e — z, u 2 )dp Ul {z) 
00 



= f P(y-0-z,u 2 )d(3 Ul (z). 



This substitution of the integral of the limit for the limit of the integral is 
legitimate because the integrands are bounded in their set. Adding and 
dividing by 2, we obtain 

3(5). i8 Ul+ u 2 (2/)= f P«,(y-z)dpuM = Pu,'Pu t (y). 

i* —00 

4. Characteristic Function. We now introduce PoincarS's characteristic 
function* which is defined as 

f " e"d0(z) = v(p, 0) Def. 

J— 00 

Before using this function it is necessary to prove some important theorems. 
* H. Poincar<§, Calcul des Probability, Paris (1912), p. 206. 
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4(1). If fi(y) is a regular function of limited variation, 

—f I [[sin k — sin (zk — yk)"2d0(z) 

exists and is equal to $(y). 
Consider the integral 

—t I [sin k — sin (zk — yk)dp(z), k 2 > ki> 0. 

The function [[sin fc — sin(zfc — yk^/irk is measurable and bounded 
(ki ^ k ^ k 2 , — oo < z < + oo ) so that the order of integration may be 
changed. If we introduce the classical symbol 

Si(x) = I dx 

Jo * 

the integral can be expressed in the form 

f^lSiik) - Si(zk - yk)JtM*)l*- 

Sinoe for real values of x, Si(x) is bounded, the integrand is limited and 
measurable, and the limit of the integral as ki = 0, k 2 = oo will exist and 
will equal the integral of the limit. (We are using again the theorem 
referred to in 3(1).) But when k = 0, Si(k) = 0, Si(zk - yk) = 0, and 
when k = oo, Si(k) = ir/2, 

Si(zk — yk) = ir/2, z> y, 
- 0, z = y, 

= - tt/2, z < y, 

so that the limit of the integrand is the idem-factor function r)(y — z). 
Therefore the double integral of the theorem exists and is equal to 

1+00 



U—u 



Corollary. — If @(y) is a regular function of limited variation and if 
<p(V> 0) = for all pure imaginary values of p = ik, then fi(y) = iden- 
tically. For then 



CO 



sin (zk — yk)d&(z) = 0, k real, 
"°dk 



^ sin fc = i/3(°o). 
Then $(y) is constant and since &(— oo) = 0, /%) = throughout. 
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4(2). If 

(** ' e h *d${z), -H<h< + K, 

J— 00 

exists, where H, K are positive, then so also does 

I z n e ht d^(z), - H < h< + K, n = 1, 2, • • • . 

*/ — 00 

For if — H < h < K we can find a number 77 > such that — H < h — 77 
< h + 57 < K. Given 77 and n we can find z > 1 such that z" < e 71 ', 

z =: z . 

Then z n e hz {z == 0) is measurable and less than the function z"«e ih+r> ' >z 
which is summable from — 00 to + °o . 

Similarly z n e hz (z ^ 0) is measurable and its modulus is less than 
z »g(A-i)s w hich is summable /8 from -co to + 00. Therefore z n e hz is 
summable /? from — 00 to + 00 . 

4(3). If 



X 



'" e h *dp(z), -H<h<K, 



exists then 



r 

t/_ or 



is holomorphic in the strip — H < real part of p < K. This strip will be 
called the HK strip. 

If p = h-\- ik lies in this strip, — H < h < K and since | cos kz | , 
I sin kz\ are not greater than 1, \e h '\ = e h ', it follows that 

X+oo /»+00 

e A * cos kz dfi(z) + i I e A * sin fcz d/3(z) 
00 J— 00; r 

exists. 

X+«o 
ze"V(2Ap)d/3(z) 

in which ^(<) = (e< - l)/<. 

Integrating along a radius I e'di = e l — 1, so that |e* — l| s |i|e'", 

and |^(0| < e |(| . 

Since — H < h < Kwe can find a number 77 > such that — H < h— 77 
< A + 77 < X and if | Ap| ^77, the real and imaginary parts of ze pz \(/(zAp) 
will not be greater in modular values than the function which has for 
each z the greater of the values of | z | e {p±1,)z , which by 4(2) is known to be 
summable. 

Furthermore yj/(t) is continuous at t = and its limit there is 1. There- 
fore by the theorem used in proving 3(1) there exists a unique limit of 
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A<p/Ap as Ap approaches 0, that is to say, <p(p) is holomorphic in the HK 
strip and 



dp 



ze" z dp(z). 

CO 



4(4). If <p(h, fa, <P(h> ft) exist, - H < h < K and if ft-ft(y) = y(y) 
then 

<p(p, 7) = <p(p, Pi)<p(P, ft) 

provided p is in the HK strip and conversely. 

<p{p, l)= f e py d v f Pi(y - z)dp 2 (z) 

«/_QO 9J—O0 

= lim e**d, piiy - z)dfo(z) 
r=» J—y «/-» 

= lim f |f e™d y pi(y - z) 1 dft(a) 

by theorem 3(1), and this 

e vz e pv dpi(v) dp 2 (z). 

x, L J- r+z J 

The expression in square brackets is not greater in modular value than 
<p(h, &>i) where h is the real part of p, «i(») is the modular variation function 
corresponding to ft(«). The existence of <p(h, Pi) implies that e hv is sum- 
mable with respect to pi, and therefore with respect to the variation function 
«i. Therefore using again the theorem with respect to limits under the 
integral sign 

<p(p, 7)= f e v Mp, Pi)dp 2 (z) 

•/ — CO 

= <p(p, Pl)<p(p, ft). 

To prove the converse let 8(y) = Pi-p 2 (y), then by the theorem itself 
<p(p, y - S) = <p{p, y) — <p{p, d) 

= <p(p, y) - <p(p, &i)<p{p, ft) 

= by hypothesis. 

Therefore by 4(1). Cor. 7 — 5=0 identically, and y{y) = /3rft(y). 

4(5). If Pu{y) is a regular function of limited variation in y for u 5: 0, 
if P{y) is not identically and if <p(p, ft) exists in the HK strip, then if 
¥>(p> ft) = e wP > where P is holomorphic in the HK strip, and where u is 
rational and non-negative, 
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ftrft,(2/) = i8«i+«,(y) 
and conversely. 

The direct theorem is a simple consequence of 4(4) since 

<p(p, A.,+0 = e (< " + " 8>p = <p( P , puMv, &,)• 

To prove the converse let p = log e <p(p, ft), i.e., where w = 1, then by 
4(4) since ft.-ft.Cy) = Pu l+ «,(y), 

log p(p, ft,) + log <o(p, ft 2 ) = log <p(p, ft 1+ u,), 

from which, if u is rational and non-negative, we obtain 

log <p(p, ft) = m log <p(p, ft) = uP, <p(p, ft) = e uP . 

To complete the theorem it remains to be proved that P is holomorphic in 
the IIK strip. Since <p(p, ft) is holomorphic P will be holomorphic except 
where <p(p, ft) = 0. Let us suppose that at p , <p(po, Pi) = then for all 
non-negative rational u, <p(po, /3 U ) = 0. But if u = nv, <p(p, j8„) = \jp(p, p v )~] n 
when n is a positive integer, and consequently S r <p(b, (5 u )2p=pi» r < n, = 0, 
where 5 stands for d/dp. 

But n can be any positive integer and r any integer less than n so that 
not only <p but all its derivatives must vanish at po. By Taylor series 
and continuation it follows that <p(p, /3„) = identically over all the IIK 
strip (boundaries excluded) and in particular that it vanishes identically 
along the imaginary axis. But by 4(1). Cor. this implies that (3 u (y) = 
identically which is a case excluded by hypothesis. This completes the 
proof of the theorem. 

This theorem enables us to pass from the functional equation 

3(5) ft.,«A., = &,+«, 

to the algebraic relation <p(p, /3 U ) = e uP and back. The formal solution 
of 3(5) can now be written in the form 

Jf* 00 dk 
-r [6" sin k - e uRik) sin { kl(k) - ky } ] 

where Rik) + il(k) = P(ik) = log <p (ik, ft), b = e m = <p(0, ft) = ft(°o) 
t« = <,«»«» =^(0, ft) = ft(*). 

5. Particular Solutions. — In our search for the more elementary solutions 
we shall carry through some processes in a purely formal manner without 
rigorous justification, but we can test the results obtained. If it is found 
that <p(p, ft) = e uP where P is holomorphic in an IIK strip then by 4(5) 
fiuiy) will be a solution of the functional equation 3(5). 
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/»00 

In some cases, particularly if I e ttS(k) dk is convergent, we may differ- 

J 

entiate both sides of 4(6) and obtain the simplified formula 

fu(y)dy 

00 

where 

f«(y) = - f'e"*** cos [*/(*) - k y yk. 

T Jo 

Let us suppose that /«(y) satisfies some differential equation which may 
be written in the form 2„ A n (D)(y n f) = 0, where D = d/dy, A n {D) is a 
polynomial in the operator D whose coefficients depend on u. Now under 
certain conditions 

f e'Vf u (z)dz = d n f e*%(z)dz, 8 = dfdp, 

= 8 n e uP . 
f " e**D[Z{z)yiz= [e**Zl±% -p f ' ' e"Zdz 

t/ — CO •-'—CO 






«/— oo 



+ 00 



C^»(D)C»y.(s!)](fe = M~ P)5 n e uP - 



Then P must satisfy the differential equation 

2U»(- P)5V P = 0. 

This equation is a differential equation for P which is a function of p only 
and therefore the ratios of the coefficients of any two distinct types of 
products of P with its derivatives must be independent of u. It follows 
that the equation can contain only two terms, those for n = and n = 1. 
Consider for example the case n = 2. Then 

M~ P)S 2 (e"0 = A 2 (- p)e nF \ju8 2 P + m 2 (5P) 2 ] 

and the ratio of the coefficients of these two terms is u and is not inde- 
pendent of u. Making the elimination of such terms there remains 

Ai(- pW + A (- p)e nP = 0, 

uA 1 (- p)5P + A (- p) = 0. 

Then^4i(— p) = A(— p) will be independent of u while A (— p) — uB(—p) 
where B is also independent of u. Since in the case considered e uPm 
= j3 M (oo) = 1, P(0) must = and 



Daniell: Integral Products and Probabilities. 155 



Pip) = £ 



A(aj d «- 



Correspondingly we have for f u (y) the equation 

A(D)(yf) + uB(D)f = 0. 

In general the simplest form of solution of this equation is given by definite 
integrals of type 4(6),* so that we do not appear to have made much 
progress. However if dP/dp is a rational function of p it may be resolved 
into a sum of terms of types a„p n , b„/(p — p ) n and if /3 M , y u correspond to 
two functions Pi, P 2 there will be a solution 5 U = j3 M -7„ corresponding to 
Pi + Pi so that more general solutions can be built out of a few simpler 
types. 

Pearsonf has suggested types of statistical distributions which are more 
general than the " normal " or Gaussian distribution. They are such 
as to satisfy a differential equation of the form 

1 df ai - y 



fdy b + biy + b 2 y 2 ' 

When this equation is cleared of fractions there appears a term b 2 y 2 Df and 
in such a case/ cannot lead to a solution of 3(5). We must therefore choose 
bi = 0, and then 

{bj) + l)(yf) + (b D - ai - h)f = 0. 
Comparing with the form already given 

A(D) = hD + 1, B(D) = (b D - |>i + b{])ju = c D - Cl . 






5(1). If bi = 0, P = icop 2 + c\P, Df/f = (ciu — y)/c u and log / 

= — y 2 /2c u + (ci/c )y + C, where C is some constant, which must be 

X+oo 
fu(y)dy = i. 
00 

In fact this gives the normal distribution 

Uy) = - 1 == «-<»-•»•***" 

■\(2ttCoU) 

in which the average is at C\U (i.e. moving with velocity c{), and the standard 
deviation is ^(cqu). 

* Cf. A. R. Forsyth, "Treatise on Differential Equations" (1903), 3d ed., p. 250. 
t K. Pearson, Philosophical Transactions of the Royal Society, 186A, p. 343. 
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<p(p, Pu) = f e™f u (y)dy 

•J— co 

1 /»+» 

= - e i<wy>Hci«P I e -(>-ClM-<;o«P) 2 /2(!oM^j. 

V(27TC M) J-*> 



= e up . 
This proves that here is a possible solution. 

5(2). If by =t= 0, let p = 1/fei, w = cipo + c pl, then 
P(p) = - c p p - n log (1 - p/po), 
log/ = — p y + (nu— 1) log (y + uc p ) + C, 

where C is some constant chosen so that /3„(oo) = 1. 

From the form of this solution it is evident that y > — uc po, or in 
other words that /„(?/) = 0, when y < — uc p . 

When y > - uc p , My) = <T P °% + wc po)" u ~ 1 p7<r u<!,,po 7r(nw). In this 
case there is a sharp boundary at y = — uc po, which is therefore moving 
with velocity — c po- The average is at (nu/p ) — uc po, the standard 
deviation is "s{nu)jpo. 

By choosing C to be partly imaginary, C = C + (raw — 1) log (— 1), 
we can obtain another solution in which f u (y) = when y > — uc po- 
For the convergence of the integral p must now be negative, = — </o 
and then when y < uc q , f u (y) = e m (uc q - y) nu - l ql u e- uc ""''lT{nu). This 
distribution is the mirror image of the other, reflection taking place at 
y = — ucoPo and with q = — po in place of po- The two cannot be com- 
bined since in one case p is positive, in the other negative. When a com- 
parison is made of the two types of solution we see that type 5(1) is con- 
tinuous and as soon as u becomes positive there is a finite non-zero density 
of distribution in both directions; type 5(2) has a discontinuity at a moving 
boundary and the distribution has non-zero density on one side only of this 
boundary. Also in 5(2) the rate of decrease in density is less rapid than 
in 5(1). To make a true comparison of these rates it is necessary to make 
allowance for the different deviations. Let y be the average y, a the 
standard deviation and put y = y + rcr where r is some pure number 
greater than 1. In case 5(1) y= C\U, <x = V(com), and 

— ffDf/f = V(c w) [cyu + r^(c u) — Cilf\/c U 
= r. 



In case 5(2) y = nu/p — uc p , a = ^(nu)/po, and 
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- <rDf/f = — - — p + 



Po L nu/p + r^(nu)/p J 

_ 1 + r V(nw) 
r + V(nw) 

If r > 1 this is < r and if r is large while u is small it is approximately I jr. 

As the next step we might consider differential equations for / which 
are of the second order, that is in which dP/dp is a fraction in which either 
the denominator or numerator (or both) is quadratic in p. If the denomi- 
nator reduces to a constant we have a case similar to 5(1) but in which 
P contains a term in p 3 . The resulting solution for / is not elementary. 
If the denominator is linear, P is the sum of terms of type 5(1) and 5(2), the 
resulting solution is an integral composition of 5(1) and 5(2) and is not 
elementary. 

Finally if the denominator is quadratic, P is the sum of types 5 (2) and 
the combination in / is again not elementary. 

There are an endless number of other possible solutions but the author 
has not had the good fortune to find types which should be so simple that 
they could be used as a basis for further investigation. Even the type 5(2), 
simple though it is, does not seem to be readily applicable to problems in 
dynamic probability. It may, however, be useful in exceptional cases. 

For the remainder of this paper we shall use the solution 5(1) in the 

forms: /3„(y) = I f u (y)dy, where 

t/ — CO 

5(3) /»(»)= - .\ e -<*-™'", 

■\(2wMu) 

5(4) f u (y) = ^(R/2wu)e- R!/y2u e Fu e- uF '' 2R , 

R = 1/M, F = RU, U = MF. 

U is the velocity of the average change in x, or we may call U the drift 
of the group. 

M is the rate of increase of the average square deviation and is called 
the mobility. The standard deviation in an interval u is \Mu. 

F is the ratio of drift to mobility and is called the force. 

R is the reciprocal of the mobility and is called the resistance. 

The formulae 5(3), 5(4) are really identical but in some cases one or the 
other form is more convenient and natural. 

6. Variable Characteristics. — We now consider the motion of a group 
during finite intervals of time and over finite intervals in the variable. 
We assume that the motion is the result of infinitesimal movements of the 
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type 5(3), that is to say, that when u is sufficiently small the motion is as 
close to that of 5(3) as we choose but that the characteristics U, M (or 
F, R) vary both with initial position x and initial time t. It is also assumed 
that, except possibly at a boundary, the first and second partial derivatives 
with respect to x of these characteristics exist and are limited, and that 
U, M, F, R are also limited. 

Let y refer to a displacement from an initial position x, u a positive 
increase in the time t and denote f„{y) by 



f(y, u; x, t) = V^ttw^-^-^ 2 ' 2 "', 
where R, U are possibly functions of x, t. Then 



i/_ o 



f(y, u; x, t)dy = 1, 



I 
£ 



f(y, u; x, t)ydy = Uu, 



f(y, u; x, t)y 2 dy = Mu, 



J J 



f(y, u; x, t)y n+1 dy = 0, n > 1, 

if we neglect terms of the second, and higher, order in u. 

In an individual case the number of members of the group lying in a 
given interval will be an integer and the distribution will be completely 
discontinuous. But our problem is that of probable, or average, not actual 
distribution and such a probable number may be fractional or even ir- 
rational. We assume that the probable number of members of the group 
in the interval x to x + dx is N(x, t)dx, that N(x, t) is limited and summable 
from — oo to + oo and that it possesses bounded first and second partial 
derivatives with respect to x. 

If K{u) is the number crossing the point x in the direction of x increasing, 

N(s, t)ds I f(y, u; s, t)dy 

= I dy I f(y, u; s, t)N(s, t)ds. 

This change of order of integration is legitimate since N(s, t) is summable 
in s and 

\f(y,u;s,t)\^^R 1 /2m- It ^e^, 
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where Ri is the maximum, R 2 the minimum value of R, Fi the maximum 
value of F = RU, the expression on the right-hand side being independent 
of s and summable in y. 

If x — p is substituted for s we obtain 

■K(«) = I «V /(J/, «; x - p, t)N(x - p, t)dp. 

Expanding in powers of p, 

f(y, u; x- p, t)N{x - p,t) = f(y, u; x, t)N(x, t) 

d p 2 

PxzU(y> u ; *> t)N(x, /)] + je(y, p), 



dx 



| e(y, p) | < max 

p=0 to y 



dx 2 



U(V> u; x- p, t)N{x - p,t)2 



< (ao + a iy 2 /u + a 2 ^/w 2 )(l + b lV + b 2 u) ^(R 1 /2u)e~ R ^ l2u e F ^, 

where ao, a\, a 2 , 61, b 2 are certain constants depending on the bounds of the 
derivatives of U, M, F, R but are independent of y and u. From this 
inequality it follows that 

I fy I nP^iy* P)dp < a term of order -w 2 . 

I tVo Jo ^ 

It can be neglected and considering terms of order u, or less, 
Z(w) = I % I /(2/> u; x, f)N(x, f)dp 

- J «"?/ J ^ I N(x, t) J f(y, u; x, t)dp 

= N(x, t) J yf(y, u; x, t)dy - ^ -^ N(x, t) J y 2 f(y, u; x, t)dy . 
Similarly the number L(u) crossing in the other direction is 
N(s, t)ds j f(y, u; s, t)dy 

x *J—co 

Xo p*-y 

dy I f(y, u; *> t)N(s, t)ds 

= - N(x, t) J yf(y, u; x, t)dy +2 fa \ N ( X > J tff(y> u > x > t) d y\ ■ 

Combining these two, the net resultant number crossing to the right (a; 
increasing) is K(u) — L(u) or 
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N(x, t) J^ yf(y, u; x, t)dy - \ §^[n(x, t) J yj(y, u; x, t)dy^ 

= NUu - ~ (NMu). 

6(1). The net rate at which members cross at x in the direction of x 
increasing is given by 

NU-H x (NM). 

6(2). Taking into consideration the flow across a neighboring point 
x + dx, if there are no sources or sinks, 

This equation is similar to that governing the flow of heat along a bar, but 
there is an added term depending on drift. 
For statistically steady flow dN/dt = 0, 

1 ^ 
NU -XT- (NM) is independent of x. 

For statistical equilibrium (statistically stationary state) 

NU-H(NM) = 

of which the general solution is N = (C/J/)e 2 * where C is an arbitrary 
constant, <p is the attractive potential defined so that d<pjdx = F = RU 
= U/M, and U, M are now assumed to be independent of t. 

It is evident that N tends to be large where the attractive potential is 
high or where the mobility is small. 

In economics <p may be regarded as a possible substitute for " desira- 
bility " or " ophelimity." 

7. Boundary Conditions. — At a rigid boundary there will not only be a 
resultant flow but the flow will be separately in each direction. Suppose 
that at x there is a distribution with finite density for smaller values of x 
and density for larger values of x; then L(u) = because N(s, t) = 0, 
s > x. Neglecting terms of order higher than u, we have seen that 

K(u) = N £ yfdy -^[n £ ffdy]. 
Let y = Uu + z a/2J/m, then 
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C yfdy = 4= P {Uu+zJ2Mrf)e-*dz 

= [ ^(Mu/2t) + \Uu~\ + terms of higher order in u. 
C X yjdy = 4= f (C/w + s <2Mufe-*dz 

Jo W J-Ut'iuIZM) 

= |Mw + terms of higher order. 
Therefore to order w, 

£(«) = l« LNU - ^ £ (iVM)H + V^VV(.af/2ir). 

This shows that the flow to the right consists of half the net statistical 
flow which would occur if the distribution were uniform through x, instead 
of being bounded at x, together with a flow in which the initial quantity 
flowing in time u is proportional to Vw. The latter constitutes a sudden 
gush from the region of finite density into the empty space. At an im- 
passable boundary this gushing flow must also vanish and the necessary 
conditions will be 

M=0, NU-~(NM) = 0. 

The condition M = makes R = <x> and this violates the conditions imposed 
in paragraph 6, so that we cannot conclude immediately that these are also 
sufficient. 

If at the boundary x , M — every member of the group which happens 
to arrive at x will have immediately thereafter a velocity U = V and the 
dispersion will be 0. It would appear to be a sufficient condition that Uo 
should be non-positive, for then the group would not pass Xq. But since 
at any point x < x 0> there is a flow of rate NU — J (d/dx)(NM) there will 
be a sudden change in the density between x and x unless this expression 
also approaches at x = x . Since M = 0, this condition is equivalent to 
Uo = \{dMjdx)o. Now M is positive except at x = x and if it possesses 
a derivative at x = x , (dM/dx) must be non-positive, or Uo "will be non- 
positive. It therefore appears that the two conditions 

M = 0, NU-~(NM)*=0 

are sufficient as well as necessary. Let W = U — ^(dM/dx), then at 
x = x , M = 0, W = 0. Now 
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At x = xo this becomes 

r) N 

~jf=- No(dW/dx) + UdM/dx) (dN/dx)o. 

It may happen that the boundary produces its effect by a " repulsion " 
which is considerable at some distance, but consider the case where U, M 
are practically unaffected until points very near x are reached. Then in 
general M will suddenly drop to from a finite value and (dM/dx) will be 
large. If (dN/dx) does not happen to be 0, dNo/dt will be large numerically 
unless there is some relation of the type dN/dx = 2kN, where & is a finite 
number such that d W/dx — JcdM/dx remains finite. Now if at the boundary 
point Xo there were statistical equilibrium without any real boundary we 
should have dN/dx = 2(W/M)N. Choose k = W'/M', where W, M' are 
the values W, M would have in the absence of a boundary. Since these 
would not vary very rapidly near Xq, the conditions would be satisfied. 
Hence we may write the general equations for our form of dynamic proba- 
bility in the form : 

dN d ,____ , 1 a 2 



7(D ~=-^ z (NU) + ^ I - 2 (NM) 



throughout the range, together with the boundary conditions : 

7(2) NU~l^(NM) = 0, 

in which U, M take the values which they would have in the absence of any 
boundary. It is important, however, to remember that in a small neighbor- 
hood of a boundary U and M actually vary rapidly and that their values 
will differ from those which they would possess in the absence of a boundary. 

Since W = = *7 - %(dM/dx) , F = R U = A^ ^r), there will 

be a strong repulsive " force " very close to a boundary. The author 

hopes to obtain some interesting results by an extension of the analysis 

to two or more dimensions. 
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